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ABSTRACT 

We present elemental abundance measurements for 9 damped Ly-a systems (DLAs) and 1 
sub-DLA at 0. 1 < z < 1.5 from recent observations with the Multiple Mirror Telescope. Most 
of these absorbers are found to be metal-poor, while 2 are found to have w 30 — 50% solar 
metallicities. Combining our data with other data from the literature, we find that the systems 
with higher [Zn/H] also have stronger depletion as measured by [Cr/Zn] and [Fe/Zn]. The re- 
lationship between the metallicity and H I column density is also investigated. Together with 
our previous MMT survey (Khare et al. 2004) we have discovered 2 of the 4 known absorbers 
at z < 1.5 that lie above (although near) the "obscuration threshold". This appears to be a 
result of selecting absorbers with strong metal lines in our sample. It would be interesting to 
find other similar systems by observing a larger sample and study how much such systems 
contribute to the cosmic budget of metals. Finally, an analysis of the iVni-weighted mean 
metallicity vs. redshift for our sample combined with data from the literature supports previ- 
ous conclusions that the iVni -weighted mean global DLA metallicity rises slowly at best and 
falls short of solar levels by a factor of > 4 even at z = 0. 
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1 INTRODUCTION 

Absorption line systems in QSO spectra provide a unique window 
into the high-redshift universe. These systems have been detected 
in the range < z < 6, corresponding to ~ 90% of the age of the 
universe. As the absorption line strengths depend primarily on the 
gas content, they are selected independent of the stellar luminosi- 
ties of the galaxies. Thus in principle they may be expected to pro- 
vide less biased probes of galaxies than objects used to study the 
high redshift universe through emission, which are biased toward 
brighter galaxies including actively star forming galaxies such as 
Lyman-break galaxies, active galaxy nuclei, and QSOs. 

Damped Lyman-a Absorption systems (DLAs) with log Nhi 
^ 20.3 and sub-DLAs with 19.0 s; log N H i < 20.3 seen in QSO 
spectra can be used as probes of the neutral gas content of the uni- 
verse. DLAs and sub-DLAs are known to contain the majority of 
the neutral gas in the universe (e.g. Wolfe et al. 1995, Peroux et 
al. 2003). With such high neutral hydrogen column densities DLAs 
are believed to be closely linked to galaxies. DLAs may also be the 



neutral gas reservoirs used in star formation. Furthermore, DLAs 
and sub-DLAs provide detailed gas-phase abundances of a number 
of elements and can thus be used to test galactic chemical evolution 
models. 

Several elements have been detected in DLAs including C, N, 
O, Mg, Si, S, Ca, Ti, Cr, Mn, Fe, Ni, and Zn. Zn is often used 
as a tracer of gas-phase metallicity as it is typically relatively unde- 
pleted in the Galactic ISM when the fraction of H in molecular form 
is low, which is the case for DLAs I Sembach et al. 1995), Zn tracks 
the abundance of Fe in Galactic stars iMishenina et al. 2002), and 
the lines of Zn II AA 2026,2062 are fairly weak and unsaturated. 
Also, Zn II can be detected with ground-based observations over 
a wide range of redshifts (0.6 < z < 3.5). Elemental abundances 
relative to Zn such as [Cr/Zn] or [Fe/Zn] can be used to estimate 
the dust content of the QSO absorber. Here and throughout the pa- 
per we adopt the standard notation [X/H] = log (Nx/Nhi) - log 
(X/H) . 

If DLAs do constitute a representative sample of galaxies, it 
would be expected that their global mean metallicity of these ob- 
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jects should rise to roughly solar values at low redshift. It is uncer- 
tain whether or not DLAs show this trend. The most recent evalua- 
tions of this by Kulkarni and Fall (2002), Prochaska et al. (2003b), 
and Kulkarni et al. (2005) show at best a weak evolution in the 
metallicity-redshift relation with slopes of ~ - 0.20±0.07. 

One proposed explanation of the lack of metal-rich DLAs in 
observed samples is that this population is being under-sampled 
due to dust obscuration (e.g.; Boisse et al. 1998; Pei & Fall 
1993; Vladilo & Peroux 2005). Absorbers with higher metallic- 
ity may also contain more dust, making any background QSOs 
fainter and keeping those objects out of magnitude-limited samples 
I York et al. 2006). In fact, it has been claimed that up to 50 % of 
QSOs may be hidden in optical surveys due to dust obscuration 
I Vladilo fc Perouxl Eofm Recently, Akerman et al. (2005) mea- 
sured the metal abundances of 2 > 1.8 DLAs from the Complete 
Optical and Radio Absorption Line System survey (CORALS), and 
found that Nhi weighted mean metallicity based on the Zn abun- 
dance is only marginally higher but within error limits of previ- 
ous samples. However, the effects of dust may be expected to be 
stronger at low redshifts. The role of dust in DLAs is still unclear, 
and more investigation is needed. We discuss dust and reddening in 
our sample of absorbers in more detail in § 5.2. 

Assuming a cosmology with £7 m = 0.3 and Qa = 0.7, z < 1.5 
spans ~70 % of the age of the universe. It is therefore critical to 
obtain more abundance measurements at low redshift. One reason 
for the uncertainty in the metallicity-redshift relation is the lack 
of abundance measurements in low-redshift absorbers. In fact ~ 
70% of the existing Zn measurements are for 2 < z < 3.5, which 
corresponds to only ~ 10% of the age of the universe. The paucity 
of Zn measurements at low z is primarily because the Zn AA 2026, 
2062 lines lie in the UV for z < 0.6, and deeply in the blue end of 
the spectrum for 0.6 < z < 1.3, where most existing spectrographs 
are inefficient. In this paper we present new MMT observations 
toward 10 QSOs. The observations presented in this paper increases 
the size of the sample of absorbers at z < 1.5 by ~ 40%. 

In § 2 we describe our observations and our data reduction 
routines. § 3 describes the estimation of column densities. In § 4 
we present notes on the individual objects. We present results from 
these measurements in § 5, and in § 6 we discuss our results and 
their implications. 



2 OBSERVATIONS AND DATA REDUCTION 

Our sample consists of 10 confirmed DLA and sub-DLA absorbers 
at 0.2 < z < 1 .5, for which Nhi is known from HST spectra 
(Rao et al. 2005). Seven of these targets were observed in the Sloan 
Digital Sky Survey (SDSS), and most show the presence of strong 
Mg II or Fe II absorption features from SDSS spectra. Apparent 
magnitudes for these 7 QSOs are available from the SDSS Quasar 
Catalog III, Data release 3 lAbazaiian et all2 0051. 

The spectra presented here were obtained at the Multiple Mir- 
ror Telescope (MMT) during two separate epochs, 2004 May and 
2005 February. Nearly 2.7 nights out of a total of 5 were lost due to 
poor weather. The Blue Chanel spectrograph was used with the 832 
1 mm - grating in the first order, or the 1200 1 mm -1 grating in the 
second order, depending on the redshift of the absorber. A CUSO4 
blocking filter was used with the 832 1 mm" 1 grating to block first 
order red light. The central wavelength for the 832 1 mm -1 grating 
was at 3670 A and 3600 A for the 2004 May and 2005 Feb epochs 
respectively. For the 1200 1 mm -1 grating the wavelength was cen- 
tered at 4935 A and 4770 A during the 2004 May and 2005 Feb 



Table 1. Atomic Data 



Species 


A-rest 

A 


/ 


References 


Bell 


3130.4219 


3.321E-1 


1 


CIV 


1548.2040 


1.899E-1 


1 




1550.7810 


9.475E-2 


1 


Mgl 


1827.9351 


2.420E-2 


1 




2026.4768 


1.130E-1 


1 




2852.9631 


1.830E00 


1 


Mgll 


2796.3543 


6.155E-1 


1 




2803.5315 


3.058E-1 


1 


A1II 


1670.7886 


1.710E00 


1 


A1IH 


1854.7184 


5.590E-1 


1 




1862.7910 


2.780E-1 


1 


Sill 


1808.0129 


2.080E-3 


1.2 


Si IV 


1393.7602 


5.130E-1 


1 




1402.7729 


2.540E-1 


1 


Can 


3933.6614 


6.267E-1 


1 




3968.4673 


3.116E-1 


1 


Till 


1910.6123 


1.040E-1 


1.3 




3383.7588 


3.580E-1 


1 


cm 


2056.2569 


1.030E-1 


1.2 




2062.2361 


7.590E-2 


1.2 


Mnll 


2576.8770 


3.610E-1 


1 




2594.4990 


2.800E-1 


1 




2606.4620 


1.980E-1 


1 


Fell 


2249.8768 


1.820E-3 


4 




2260.7805 


2.440E-3 


4 




2344.2139 


1.140E-1 


1.5 




2374.4612 


3.130E-2 


1 




2382.7652 


3.200E-1 


1.5 


Znll 


2026.1370 


5.010E-1 


1,6 




2062.6604 


2.460E-1 


1.6 



References. (1) Morton 2003; (2) Bergeson & Lawler 1993b; 

(3) Wiesse et al. 1996; (4) Bergeson et al. 1994; (5) Bergeson et al. 1996; 

(6) Bergeson & Lawler 1993a 



epochs respectively. In order to achieve the desired S/N while min- 
imizing cosmic rays, multiple exposures of each target were taken 
with exposure times ranging from 1800 to 2700 s depending on 
the magnitude of the QSO. Each target observation was followed 
by a comparison spectrum from a He+Ne+Ar lamp for wavelength 
calibration. Quartz flat fields and bias frames were taken at both 
the beginning and end of the night. Table 1 lists the vacuum wave- 
lengths and oscillator strengths used in identification of the features 
and subsequent analysis.Table 2 gives a summary of our observa- 
tions. 

Data reduction was carried out using standard IRAF routines. 
Images were bias subtracted, trimmed, and flat fielded using the 
task CCDPROC. These reduced 2-d spectra were then extracted 
and wavelength calibrated using the task DOSLIT The wave- 
lengths were converted to vacuum using the task DISPTRANS. Fi- 
nally the spectrum was normalized using the task CONINUUM. 
The solar abundances used as reference while calculating the abun- 
dances are adopted from Lodders et al. (2003). 

Figures 1-10 plot the co-added and continuum-fitted spectra 
obtained with the MMT, along with line identifications from the 
DLAs as well as other absorption systems along the line of sight. 
Unmarked, minor depressions have been verified to be artifacts. 
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Table 2. Summary of Observations 



QSO 


9 


Zern 


ZDLA 


log N H i 

_2 

cm 


Grating 
lines mm 


FWHM 

A 


Exposure Time 


Epoch 


0738+313 


16.1 


0.630 


0.0912 


21.18±0.06 


832 


1.10 


3000 


2005 Feb. 








0.2210 


20.90±0.09 


832 


1.10 


3000 


2005 Feb. 


0827+243 


17.3 


0.941 


0.5247 


20.30±0.05 


832 


1.10 


3600 


2005 Feb. 


1010+0003 


18.3 


1.398 


1.2651 


21.52±0.07 


1200 


1.33 


3600 


2005 Feb. 


1107+0003 


18.7 


1.741 


0.9547 


20.26±0.14 


832 


1.10 


8100 


2005 Feb. 


1137+3907 


17.4 


1.023 


0.7190 


21.10±0.10 


832 


1.10 


5400 


2004 May 


1225+0035 


18.9 


1.226 


0.7731 


21.38±0.11 


832 


1.10 


5400 


2005 Feb. 


1501+0019 


18.1 


1.930 


1.4832 


20.85±0.13 


1200 


1.33 


5400 


2004 May 


1712+5559 


18.7 


1.356 


1.2095 


20.72±0.05 


1200 


1.33 


5400 


2004 May 


1715+4606 


18.1 


0.989 


0.6511 


20.44±0.13 


832 


1.10 


5400 


2004 May 


1733+5533 


18.1 


1.074 


0.9981 


20.70±0.03 


832 


1.10 


3600 


2004 May 



3 ESTIMATION OF COLUMN DENSITIES 

Column densities were estimated by fitting the observed absorption 
line profiles using the package FITS6P l W eltv et allll99lh . which 
has evolved from the code by Vidal-Madjar et al (1977). FITS6P 
iteratively minimizes the % 2 value between the data and a theoreti- 
cal Voigt profile that is convolved with the instrumental profile. All 
lines were fit with a single component Voigt profile. If a multiplet 
was available such as Fe II AA 2344, 2374, 2382 the lines were 
fit simultaneously until convergence for a common N, 6 e //, and 
V. Equivalent width measurements were obtained with the package 
SPECP, also written by D.E. Welty. The la equivalent width un- 
certainties were estimated using the photon noise uncertainties and 
the continuum uncertainties obtained by allowing the continuum to 
vary by ± 10 %. In all cases these uncertainties were dominated by 
the photon noise uncertainties. 

The prescription of Khare et al. (2004) for fitting the Cr II and 
Zn II lines of was followed. Specifically, the Cr II AA2056, 2066 
lines were fit simultaneously to obtain N, 6 e //, and v values. Then, 
the blended Cr Il+Zn II A 2062 line was fit starting with the values 
obtained from the previous fits, holding the Cr II component fixed. 
The Zn II A 2026 line was then fit using the b e ff and v value from 
the Cr AA 2056, 2062 fit while holding the Mg I component fixed. 
If the Cr II A 2056 line could not be fit due to noise, we fit the Zn 
II A2026 line, holding the Mg I A 2026 contribution fixed, then fit 
the Cr II + Zn II A 2062 blend holding the Zn II component fixed. 
The Mg I column density was estimated by fitting the Mg I A 2852 
line from the SDSS spectra. York et al. (2006) showed that the Mg 
I A 2852 lies on the linear portion of the curve of growth for W^gi 
~ 0.6 to 0.7. The Mg I A 2852 line is often moderately saturated 
in DLAs and the assumption of a linear curve of growth may lead 
to an underestimation of the Mg I column density, therefor we did 
not assume Mg I A 2852 to be on the linear portion of the curve 
of growth, and instead used the results of a profile fit to that line. 
To judge the effects of possible saturation of Mg I A 2852 on NMgi 
and hence on Nznii, we also estimated the maximum NMgi that 
would be consistent with W2852 within a 3<r level. We then used 
this maximum NMgi to estimate the minimum Nznii in the A 2026 
line. In any case, in our experience the contribution of Mg I to the 
blended A 2026 line is typically small (eg. Peroux et al. 2006a). 
For Q 150 1+00 19 we also have a constraint on NMgi from the non- 
detection of the Mg I A 1827 line which gives log NMgi < 13.47. 

Tables 3 and 4 list the derived equivalent widths and column 
densities respectively. If a certain species was not detected, 3<r up- 
per limits were placed on the equivalent widths and column densi- 



ties based on the signal to noise ratio (S/N) in the region based on 
a 3 pixel resolution element. It was assumed that these lines lie on 
the linear portion of the curve of growth. As can be seen in figures 
1-10, none of the detected Zn II or Cr II lines are saturated. In fact, 
none of the lines detected appear to show significant saturation even 
in those cases where the equivalent widths are large. This is likely 
due to the line profile being composed of multiple components that 
could not be resolved at the resolution of our data. 



4 DISCUSSION OF INDIVIDUAL OBJECTS 

Q0738 + 313 (z em = 0.630, z abs = 0.0912 for system A and z abs 
= 0.2210 for system B): Both of the systems listed are DLAs. The 
Zn and Cr lines were studied by Kulkarni et al. (2005) using HST 
data. The AA 3933,3969 lines of Ca II appear in both systems A 
and B. We also placed upper limits on Ti II A 3383 of system B. 
As the background QSO is very bright (g = 16.1), we were able to 
obtain high S/N of ~ 100 in the region. This gives an upper limit of 
log (Ntiii/Nhi) < -9.42, with [Ti/H] < -2.32. The feature located 
at ~ 5070 A is an artifact from bad pixels in the CCD. 

Q0827 + 243 (z em = 0.941, z abs = 0.259 for system A and 
z abs = 0.5247 for system B): The Zn and Cr lines of system B 
were studied in detail at lower wavelengths by Khare at al. (2004). 
The Fe II AA 2249, 2260, 2344, 2374, 2382, 2586, 2600 lines were 
measured previously. Strong Mg II AA 2796, 2803 and Mg I A 2852 
lines were present in system B in the new spectra. The Mg I A 2852 
line was fit first to estimate b e // and v values, which were used 
when fitting the Mg II doublet. Due to the strong saturation of the 
Mg II AA 2796, 2803 doublet we give the derived column density 
as a lower limit. We were able to place upper limits on Ti II A 3383 
of [Ti/H] = - 1 .45 at S/N ~ 60 in the region. The feature located at 
~ 5070 A is again an artifact from bad pixels in the CCD. 

Q1010 + 0003 (Zem = 1.4007, z abs =1.2651): The Cr II AA 
2056, 2062, 2066 lines, as well as the Zn II AA 2026, 2062 lines 
were detected with S/N of ~ 30. As all the lines were detected, the 
full prescription outlined in § 2 was followed. The Mg I contribu- 
tion to the blended Zn II+Mg I A 2026 line was estimated from the 
SDSS spectrum. The Mg I A 2852 has an equivalent width of W res t 
= 404 mA. This line was fit using a single component model, yield- 
ing log NMgi = 12.67±0.05. This component was held fixed in Mg 
I A 2026 while the Zn II and Cr II AA 2026, 2056, 2062, 2066 lines 
were fit simultaneously giving log Nznii = 12.96. 
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Table 3. Equivalent Widths 



QSO 


'' a b s 


Species 


Arest 

A 


Wrest 

mA 


QSO 


Zabs 


Species 


^rest 

A 


Wrest 

mA 


0738+313 


0.0912:A 


Call 


3933 


189±13 


1501+0019 


1.4832 A 


Mgl 


2852 


952 a 


Zem = 0.630 




Call 


3969 


54±10 


z em = 1.930 




Aim 


1854 


334±17 




0.2210:B 


Call 


3933 


63±10 






Aim 


1862 


216=1=18 






Call 


3969 


35±10 






Sill 


1808 


295±16 






Till 


3383 


<14 






CrII 


2056 


105±16 


0827+243 


0.2590:A 


Call 


3933 


<34 






CrII 


2066 


<26 


z em = 0.941 




Call 


3969 


<34 






Zn II+Mg I 


2026 


231=1=21 






Tin 


3383 


<33 






Zn Il+Cr U 


2062 


191=1=18 




0.5247:B 


Mgl 


2852 


602±29 




1.8510:B 


CIV 


1548 


154=1=18 






Mgll 


2796 


2416±27 






CIV 


1550 


116=1=18 






Mgll 


2803 


2316±28 




1.9260:C 


CIV 


1548 


254=1=14 






Till 


3383 


<30 






CIV 


1550 


154=1=14 


1010+0003 


1.2651:A 


Mgl 


2852 


404 a 






Ain 


1670 


88=1=11 


Zem = 1-398 




A1III 


1854 


104±21 


1712+5559 


1.1590:A 


CrII 


2056 


<33 






Aim 


1862 


<46 


Zem == J--OOD 




CrII 


2066 


<33 






CrII 


2056 


91±19 






Fell 


2344 


415±30 






CrII 


2066 


69±15 






Fell 


2374 


240±27 






Fell 


2249 


144±26 






Fell 


2382 


592±32 






Fell 


2260 


184±26 






Zn II+Mg I 


2026 


<33 






Fell 


2344 


564±28 






Zn Il+Cr II 


2062 


<33 






Fell 


2374 


543±35 




1.2093:B 


Mgl 


2852 


367° 






Fell 


2382 


715±30 






CrII 


2056 


<35 






Zn II+Mg I 


2026 


233±23 






CrII 


2066 


<35 






Zn Il+Cr II 


2062 


181±22 






Fell 


2344 


973±33 


1107+0003 


0.5252:A 


Mnll 


2576 


<30 






Fell 


2374 


723 ±31 


Zem = 1.741 




Mnll 


2594 


<34 






Fell 


2382 


1168=1=35 






Mnll 


2606 


<20 






Zn II+Mg I 


2026 


<35 






Fell 


2344 


139±36 






Zn Il+Cr H 


2062 


<35 






Fell 


2374 


<47 


1715+4606 


GAL 6 


Call 


3933 


1556±46 c 






Fell 


2382 


269±36 


Zem = 0.989 




Call 


3969 


205±36 






Fell 


2586 


185±25 




0.6544A 


CrII 


2056 


<165 






Fell 


2600 


342±27 






CrII 


2066 


<165 




0.9547:B 


Tin 


1910 


<42 






Fell 


2344 


1264±42 






CrII 


2056 


<27 






Fell 


2374 


851 ±40 






CrII 


2066 


<27 






Fell 


2382 


1556±46 c 






Zn II+Mg I 


2026 


<27 






Zn II+Mg I 


2026 


<165 






Zn Il+Cr H 


2062 


<27 






Zn Il+Cr H 


2062 


<165 




1.711:C 


Si IV 


1393 


618±36 


1733+5533 


GAL 6 


Call 


3933 


482±36 






Si IV 


1402 


447±40 


z em = 1.074 




Call 


3969 


257±33 


1137+3907 


GAL 6 


Call 


3933 


336±58 






Mgl 


2852 


362 a 


Zem — — l.Uzo 




Call 


3969 


241±51 




0.9984A 


Aim 


1854 


179=1=31 




0.7190:A 


CrII 


2056 


<74 






Aim 


1862 


100±27 






CrII 


2066 


<74 






Sill 


1808 


178±30 






Fell 


2249 


150±31 






Till 


1910 


<30 






Fell 


2260 


291±34 






CrII 


2056 


<35 






Fell 


2344 


1768±67 






CrII 


2066 


<29 






Fell 


2374 


1300±37 






Zn II+Mg I 


2026 


<29 






Fell 


2382 


2009±70 






Zn Il+Cr H 


2062 


<29 






Zn II+Mg I 


2026 


447±56 




1.1496:B 


CIV 


1548 


453±55 






Zn Il+Cr H 


2062 


266±70 






CIV 


1550 


341=1=51 


1225+0035 


0.7731 A 


Mgl 


2852 


929" 






Al II 


1670 


<40 


Zem — — i-.£Z\y 




CrII 


2056 


259±81 






Aim 


1854 


<30 






CrII 


2066 


177±76 






Aim 


1862 


<30 






Fell 


2249 


339±73 






Sill 


1808 


<25 






Fell 


2260 


372±80 
















Zn II+Mg I 


2026 


<122 
















Zn Il+Cr H 


2062 


181±91 













a : From SDSS spectra. b : Entries with GAL identifier are from lines originating in the Milky Way. c : The Galactic Ca II A 
3933 line and Fe II A 2383 line from the DLA at z a j, s = 0.6544 are blended in the spectra. 
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Figure 1. Combined and normalized spectrum of Q0738+313. The expected positions of several lines of interest are marked above the continuum with their 
rest wavelength and species. The letters after the wavelengths refer to the different systems indicated at the top of the figure. An "X" signifies a defect in the 
CCD. 
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Figure 2. Same as Figure 1, but for Q0827+243. 
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Figure 8. Same as Figure 1, but for Q1712+5559. 
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Figure 9. Same as Figure 1, but for Q1715+4606. 
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Figure 10. Same as Figure 1, but for Q1733+5533. 
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To judge the effect of potential saturation of Mg I A 2852 we 
estimated the maximum NMgi and hence the minimum Nznii as 
described in § 3. This gave log Nznii > 12.95, affirming the as- 
sumption that the Mg I contribution to the blended A 2026 line is 
small. Fe II AA 2249, 2260, 2344, 2374, 2382 lines were also de- 
tected. 

Q1107 + 0003 Oem = 1.7408, z abs = 0.5252 for system A, 
Zabs = 0.9547 for system B and z abs = 1.71 1 for system C): System 
B is of primary interest to this paper, and contains a confirmed sub- 
DLA from HST spectra with log Nhi = 20.26. This QSO also has a 
broad absorption line system (BAL) at z ~ 1.711. Si IV AA 1393, 
1403 lines from the BAL were detected. No Zn II or Cr II lines 
were detected from the sub-DLA at z = 0.9547, but upper limits 
were placed on the column densities. There is also a system at z 
= 0.5252 with weak Fe II lines detected in SDSS spectra. Lines of 
Fe II AA 2344, 2382, 2586, 2600 were detected from this system. 
The weaker Fe II A 2374 line however was not detected. The Fe 
II A A 2344, 2382 lines were fit simultaneously, as were the AA 
2586, 2600 lines with the same beff and v. Mn II AA 2576, 2594, 
2606 fell within the observable region for system B, but none were 
detected. 

Q 1 137+3907 (z em = 1.023 and z abs =0.7193): This is aBAL 
system. The Zn II + Mg I A 2026 blend was detected at ~ 3cr level 
in this system. The S/N in the region was ~ 15, leading to larger 
errors in the column density estimates. There are no data available 
on Mg I for this system, so it could not be included in the fit of the 
blended Zn II+Mg I A 2026. We are however confident of the Zn II 
column density quoted, because based on our experience so far, the 
contribution from Mg I A 2026 is relatively small for most systems. 
For example, the z = 1.2651 DLA toward Q 10 10+0003 showed 
differences in the Zn II column density of only ~ 0.05 dex with 
and without the contribution from Mg I A 2026.The Zn II A 2026 
line was fit using the beff value from the fit of the Fe II AA 2249, 
2260 lines. The Zn II component was held fixed in the blended A 
2062 line to determine Ncni. Strong Fe II A A 2249, 2260, 2344, 
2374, 2382 were also detected. The feature at ~ 3740 A is caused 
by a defect in the CCD. 

Q1225 + 0035 (z em = 1.226 and z abs = 0.7730): The spectra 
of this object were taken near daybreak, which resulted in addi- 
tional noise due to increased background. The Cr II AA 2056, 2062 
lines was detected at ~ 3a, allowing for moderately accurate deter- 
mination of the column density. The Zn II+Mg I A 2026 blend lies 
in a noisy region and was not detected, and the Zn II upper limit 
had to be estimated from S/N. We also detected Fe II A 2249, 2260 
lines which were fit simultaneously with the same beff and v. The 
feature at ~ 3735 A is caused by a defect in the CCD. 

Q1501 + 0019 (z em = 1.930, z abs = 1 .4832 for system A, z abs 
= 1.8510 for system B, and z abs = 1.9260 for system C): System A 
is the system of primary interest in this paper, and has a confirmed 
DLA from HST spectra. We detected only the C IV A A 1548, 1550 
lines in system B. The Zn II+Mg I A 2026 blend, along with the Cr 
Il+Zn II A 2062 and Cr II A 2056 lines were all detected in system 
A. Because the Cr II A 2066 line was not detected, the Cr II column 
density was determined from Cr II A 2056 solely. The Mg I contri- 
bution to the blended Zn II+Mg I A 2026 line was determined from 
SDSS spectra, which yields a W rest for the Mg I A 2852 line of 
952 mA. This line was fit using a one component model, giving log 
N Mg i = 12.93±0.05. This component was held fixed in the MMT 
spectrum of the A 2026 line while the Zn II and Cr II AA 2056, 
2062, and 2066 lines were fit simultaneously. To judge the effects 
of saturation of Mg I A 2852, we again estimated the maximum Mg 



I contribution to the blended A 2026 line as described in § 3. This 
gave Nznii > 12.89, only 0.04 dex less than the "best fit" model. 
Al III A A 1854, 1862 lines were also detected and fit together while 
simultaneously varying N, b e ff, and v. The feature at ~ 5220 A is 
caused by a defect in the CCD. 

Q1712 + 5559 (z em = 1.930, z abs = 1.2093 for system A 
and z a bs = 1.1584 system B): System A is the system of primary 
interest for this paper, and contains a confirmed DLA from HST 
spectra with log Nhi = 20.72. Both systems contain strong Fe II 
AA 2344, 2374, 2382 lines. No other lines from system B appear to 
be present. No Zn II or Cr II lines were detected in either system, 
with S/N ~ 30 in the region. The feature at ~ 5220 A is caused by 
a defect in the CCD. 

Q1715 + 4606 (z em = 0.989 and z abs = 0.6511) We have 
determined the H I column density for this object from archival 
HST data. By fitting the line profile, we have determined a column 
density of log N H i=20.44. Strong Fe II AA 2344, 2374, 2382 lines 
were detected. No lines of Zn II or Cr II were seen, although the 
S/N in the region was only ~ 10 due to poor weather during the 
observations. A weak Galactic Ca II A 3698 feature was seen, but 
the stronger Ca II A 3933 line was blended with the Fe II A 2382 
line from the z abs =0.65l 1 system. The feature located at ~ 3690 A 
is an artifact of the CCD. 

Q1733 + 5533 (z em = 1.074 and z abs = 0.9984 for system A 
and Zabs = 1.1496 for systemB) System Ais of primary interest for 
this study, as it has a confirmed DLA with log Nhi = 20.70. System 
B appears to be a C IV system with no other absorption features 
at this redshift detected in SDSS spectra. No Zn II or Cr II lines 
were detected in system A, but the relatively high S/N (~33) in the 
region allows us to place tight upper limits on the abundances. Al 
III A A 1854,1862 lines were weakly detected and were fit together 
varying N, 6 e //, and v simultaneously. We also detected strong 
Galactic Ca II A A 3933, 3968 lines. 



5 RESULTS AND DISCUSSION 

5.1 Relative Abundances 

The observed relative abundances of the elements are a combina- 
tion of both the nucleosynthetic processes and of dust depletion. 
If dust is significantly present in DLAs, then refractory elements 
such as Cr, Mn, Fe, Co, and Ni should show substantial depletions. 
As already discussed in § 1, Zn is a relatively undepleted element 
due to its low condensation temperature. Ratios such as [Cr/Zn] 
and [Fe/Zn] are therefore a measure of dust depletion. Table 5 lists 
abundances of Cr, Fe, and Ti relative to Zn for the systems in our 
sample. Also given is the metallicity [Zn/H]. Although it is cus- 
tomary to estimate metallicity based on measurements of Nznii, 
it has b een shown that [Zn/H] is not much affected by ionization 
IVladilo et alfeOOll) . 

No Ti II was found in any of the systems sampled. We deter- 
mined [Ti/Zn] upper limits of -1.01 and -1.24 for the DLAs toward 
Q1010+0003 and Q1501+0019 respectively. The measured [Fe/Zn] 
in Q1010+0003 was also sub-solar (-0.69) also suggesting dust de- 
pletion. No Fe lines were sampled in the spectra of Q1501+0019, 
so [Fe/Zn] could not be measured. 

Si is an a-capture element, and its abundance relative to Fe 
group elements can provide clues into the chemical evolution of 
DLAs. Disk stars in the Milky Way show a systematic decrease 
in [a/Fe] with increasing [Fe/H], which indicates an increase with 
time of Fe contributed to the ISM f rom type la supe rnovae, rel- 
ative to type II supernovae JEdvardsson et alJ ll993V Prochaska 
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Table 4. Column Densities 
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15.71±.02 


86.5 






Till 


<11.48 








CrII 


13.40±0.09 


86.5 






CrII 


13.11±0.24 6 








Znll 


12.93±0.06 


86.5 






Znll 


<12.83 b 






1.8510 


CIV 


13.64±0.04 


74.1 


0827+243 


0.2590 


Call 


<11.60 






1.9260 


CIV 


13.87±0.03 


88.8 






Till 


<11.96 








A1II 


12.33±0.04 


53.4 




0.5247 


Mgl 


12.69±0.02 


80.3 


1712+5559 


1.1584 


CrII 


< 12.93 








Mgll 


>14.59 


80.3 






Fell 


13.98±0.03 


74.9 






Till 


<11.76 








Znll 


< 12.25 








CrII 


<13.42 6 






1.2093 


Mgl 


12.44 c 


76.7 






Fell 


14.59±0.02 a 








CrII 


<12.86 








Znll 


<12.80 6 








Fell 


> 14.54 


76.7 


1010+0003 


1.2651 


Mgl 


12.67±0.05 c 


184.3 






Znll 


<12.19 








A1III 


12.67±0.10 


10.8 


1715+4606 


GAL d 


Call 


12.66±0.10 


35.5 






CrII 


13.54±0.07 


22.5 




0.6511 


CrII 


<13.54 








Fell 


15.26±0.05 


22.5 






Fell 


14.94±0.03 


57.4 






Znll 


12.96±0.06 


22.5 






Znll 


<12.87 




1 107+0003 


0.5252 


Mnll 


<12.16 




1733+5533 


GAL d 


Call 


12.79±0.03 


88.2 






Fell 


13.54±0.08 


37 




0.9984 


Mgl 


12.44 b 






0.9547 


Till 


<13.01 








A1III 


13.13±0.05 


133.7 






CrII 


<12.76 








Sill 


15.48±0.06 


118.9 






Znll 


< 12.08 








CrII 


< 12.79 






1.7110 


Si IV 


13.99±0.03 


198.3 






Znll 


<12.11 




1137+3907 


GAL d 


Call 


12.63±0.05 


87.8 




1.1496 


CIV 


14.19±0.04 


88.8 




0.7193 


CrII 


13.71±0.20 


122.9 






A1II 


<11.77 








Fell 


15.45±0.05 


122.9 






SiH 


<14.6 








Znll 


13.43±0.05 


122.9 













a : Khareet al. (2004) b : Kulkarni et al. (2005) c : FromSDSS spectra. d : Entries with GAL identifier are from lines originating 
in the Milky Way. 



& Wolfe (1999) observed that the DLAs in their sample showed 
no such trend, and that most of the objects were a-enhanced. In 
DLAs, Zn is a better indicator of Fe group abundances than Fe 
because of the large depletion of Fe on dust grains. We detected 
Si II A 1808 in Q1501+0019 and Q1733+5533. [Si/H] = -0.68 in 
Q1501+0019, and [Si/H] = -0.76 in Q1733+5533. We unfortunately 
did not cover any Fe II lines in the region. However, [Si/Zn] = -0.28 
in Q1501+0019, and [Si/Zn] > 0.37 in Q1733+5533, implying sig- 
nificant a-enhancement in the latter system. 

Figure 11 shows [Cr/Zn] vs. [Zn/H] for our data as well as 
those from the literature. Only detections, not limits have been 
plotted. There is intrinsic scatter in the data, which is to be ex- 
pected, but there appears to be a general trend suggesting that the 
most metal rich DLAs also show the most depletion. If this is 
true, then one would expect the most metal rich absorbers to be 
under-sampled because they would also most likely be the most 
reddened. A Spearman rank-order correlation test on the 48 data 
points plotted in Fig. 1 1 gives the Spearman rank-order coefficient 
R s =-0.623, with a probability of obtaining such a value by chance 
ofP(R s )=1.95xlO~ 5 . 

Figure 12 shows [Fe/Zn] vs [Zn/H] for our data as well as 
from the literature. Again, only detections have been plotted. Both 



figures 11 and 12 show a trend toward higher depletion with 
higher metallicity. Previous investigations have shown a similar 
anti-correlation between [Cr/Zn] or [Fe/Zn] and [Zn/H] although 
the vast majority of the i r data set has come from absorbers at z > 
1.5 jPettini et alJll997l : IProchaska fe Wolfell2002l lAkerman et all 
2005). Our study increases the number of Zn measurements at 
z < 1.0 by ~ 30% A Spearman rank-order correlation test on 
this data set gives Rs=-0.428, with a probability of obtaining such 
a value by chance of P(Rs)=6.00x 10~ 3 . 



5.2 Dust Obscuration Bias 

Several studies have tried to understand whether current DLA ob- 
servations are affected by a dust obscuration bias. Recently, York et 
al. (2006) conducted a survey of 809 Mg II systems from the SDSS 
with 1.0 ^ z ^ 1.9 and found that QSOs with an intervening ab- 
sorber are at least three times more likely to have highly reddened 
spectra than QSOs without any absorption systems in their spec- 
tra. Furthermore, York et al. (2006) found a trend toward higher 
Zn II column densities but lower [Zn/H] and lower abundances of 
the first ions as the reddening E(B-V) increased. Significant red- 



Table 5. Relative Abundances 
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QSO 


^abs 


[Zn/H] 


[Cr/Zn] 


[Fe/Zn] 


[Ti/Zn] 


[Si/Zn] 


Q0738+313 


0.0912 


<-1.14 a 


>-0.21 a 


>-0.48 a 








0.2213 


<-0.7 a 


<0.7 a 










Q0827+243 


0.5247 


<-0.04 a 


<-0.3 a 










Q1010+0003 


1.2651 


-1.19 


-0.46 


-0.54 


<-1.01 






Ql 107+0003 


0.9547 


<-0.72 












Ql 137+3907 


0.7190 


-0.30 


-0.67 


-0.82 








Q1225+0035 


0.7731 


<-1.15 


>0.09 


>0.07 








Q1501+0019 


1.4832 


-0.54 


-0.56 




<-1.24 


-0.28 


Q17 12+5559 


1.2093 


<-1.07 




>-0.51 






Q1715+4606 


0.6511 


<-0.61 




>-0.78 






Q1733+5533 


0.9981 


<-1.13 








>0.37 



1 : Abundance measurements from Kulkarni et al. (2005). 
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Figure 11. [Cr/Zn] vs. [Zn/H] for our data as well as measurements from 
the literature. Only detections are plotted. 
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Figure 12. [Fe/Zn] vs. [Zn/H] for our data as well as from the literature. 
Only detections are plotted. 



dening has also been seen in a s ample of strong Ca II absorbers 
IWild fe Hewetj2005tlwiid et all200d) . 

If dust obscuration effects are significant, highly reddened 
QSOs would be under-sampled due to their faintness. Boisse et al. 
(1998) pointed out that there is a deficit of systems with high iVm 
and high [Zn/H], and suggested that this may be caused by a dust 
obscuration bias (because such systems may also have high dust 
content and may obscure the background quasars more). They spec- 
ulated that observations of fainter QSOs might reveal systems with 
high TVhi and high metallicity. Figure 13 shows a plot of [Zn/H] vs. 
iVm for DLAs from our sample as well as the literature (Akerman 
et al. 2005, Rao et al. 2005). Clearly, there is a deficit of systems 
that are seen with both low metallicity and low 7Vm. This could be 
attributed to observational limitations in detecting weak Zn II lines, 
although recent high resolution, high S/N observations still find a 
dearth of systems in this region. More interestingly though, there 
are few absorbers with high metallicity and high Nnx, as noted by 
Boisse et al. (1998) and also by Akerman et al. (2005). This can- 
not be credited to observational limitations, because systems with 
similar metallicities have been detected at lower Nja, an d seems to 
suggest that these types of systems are being under-sampled. Also, 
there is nothing unphysical about systems with high TVhi and high 
metallicity. It is thus surprising that only 3 points from the literature 
(of which 1 came from our previous sample in Khare et al. 2004) lie 
above the dashed line indicating the empirical "obscuration thresh- 
old" log iVznii < 13.15 suggested by Boisse et al. (1998). 

Two of the systems at z < 1.5 that lie above the "obscura- 
tion threshold" came from our samples (this paper and Khare et al. 
2004), although it is still remarkable that all points in Fig. 13 lie 
fairly close to the threshold suggested by Boisse et al. (1998). In- 
deed, all but two of the points are consistent, within the errors, with 
being below the obscuration threshold. It should also be noted that 
the Boisse et al. threshold is not meant to be a hard limit, but an 
observationally determined boundary based on the much smaller 
sample available to Boisse et al. (1998). Indeed, hydrodynamical 
simulations by Cen et al. (2003) predict the existence of absorbers 
above this threshold. Thus, it may not be completely surprising that 
we are now finding some objects above the threshold. 

We emphasize that the higher proportion of systems with 
Nznii close to the Boisse et al. threshold in our studies (this paper 
and Khare et al. 2004) is not due to errors in the column density es- 
timations from the moderate resolution of the MMT spectrograph. 
Column densities derived from our observations of several systems 
with both the MMT and with the higher resolution (~ 5 km s _1 ) 
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Ultraviolet- Visual Echelle Spectrograph (UVES) on the Very Large 
Telescope (VLT) agree to within 0. 1 dex JPeroux et all2006ala) . A 
similar agreement was also found between column densities de- 
rived by Khare et al. (2004) from MMT spectra and those derived 
by Rao et al. (2005b) from Keck HIRES spectra. The reason for the 
higher proportion of metal-rich systems in our studies is that the 
targets in this paper and Khare et al. (2004) were chosen partially 
because of strong metal line absorption features seen in SDSS spec- 
tra. Thus, these systems may have been more likely to have higher 
Zn column densities. Indeed, composite SDSS quasar spectra indi- 
cate that systems with larger W res t of Mg II A2796 tend to have 
stronger Zn II lines (Nestor et al. 2003; Yorketal. 2006). One of our 
two systems has a large depletion consistent with the obscuration 
selection effect while the other system shows a moderate depletion, 
so it is not clear whether the small number of such systems found 
so far is caused by dust obscuration or small sample size especially 
at z < 1.5 where the effects of dust are expected to be the most 
pronounced (Boi sse et alJ ll998). It could also be the case that the 
most metal rich systems are truly rare. 

Figures 14 and 15 show the dust content as measured by 
[Cr/Zn] and [Fe/Zn] respectively vs. log Nm for these DLAs as 
well as those from the literature. As can be seen, there is very little 
trend in the data. It may appear surprising from Figs. 11-15 that 
there is no (anti)-correlation between [Cr/Zn] (or [Fe/Zn]) and log 
N(HI) even though there is an (anti)-correlation between [Cr/Zn] 
(or [Fe/Zn]) and [Zn/H] and between [Zn/H] and log Nhi- If one 
traces some of the extreme points in Figs. 1 1 (or 12) and 13 that do 
show the correlation, one can see that the correlations in these two 
figures arise from different systems. For example, consider the 3 
points A, B, and C with ([Cr/Zn], [Zn/H], log Nhi « (-0.46, -1.19, 
21.52); (-0.74, -0.30, 21.10); and (-0.54; -.57; 20.85). Points such 
as A and B are responsible for the correlation in Fig 1 1 . Points such 
as B and C are responsible for the correlation in Fig. 13. When one 
plots these points in Fig. 14, there is no strong correlation. (More 
generally, if data sets x,y and y,z show correlations, x and z do 
not necessarily show a correlation; whether or not they do depends 
on the distribution of the individual data values. See, for example, 
Casella & Berger 2002.) 



5.3 Constraints on Metallicities and Potential Implications 
for Metallicity Evolution 

The Zn abundances for our DLAs are listed in Table 5. In 4 of the 
systems at 0.6 < z < 1.5 the abundances of Zn are about 10% 
solar or lower, while in 2 systems, the abundances are 30-50 % 
solar. We now briefly discuss the implications of our data for the 
metallicity evolution of DLAs. Our analysis is based on the statis- 
tical procedures outlined in Kulkarni & Fall (2002), and uses our 
data, the data compiled in Kulkarni et al. (2005), and other recent 
data from the literature (Rao et al. 2005; Akerman et al. 2005). We 
binned the combined sample of 109 DLAs in 6 redshift bins with 18 
or 19 systems each and calculated the global Nhi -weighted metal- 
licity in each bin. Figure 16 shows [Zn/H] vs. z for the data in the 
literature as well as from this paper. 

We constructed two samples for these 109 DLAs: a "maxi- 
mum limits" sample where the Zn limits are treated as detections, 
and a "minimum limits" sample, where the Zn limits are treated 
as zeros. For an individual system these extreme cases cover the 
full range of possible values the Zn column densities can take in 
the case of the limits. The iVm-weighted mean metallicity in the 
lowest redshift bin 0.1 < z < 1.2 is -0.78 ± 0.11 for the "max- 
imum limits" sample and —0.99 ± 0.17 for the "minimum limits" 
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Figure 13. [Zn/H] vs. log Nuj for our data as well as previous measure- 
ments from the literature. The dashed line indicates the empirical upper 
limit log -/Vznll < 13.15 inferred from previous studies and suspected to 
be the "obscuration threshold". Only detections are plotted. 



sample. The linear regression slope of the metallicity-redshift re- 
lation for the redshift range 0.1 < z < 3.9 is -0.18 ± 0.07 for 
the "maximum limits" sample, and —0.25 ± 0.10 for the "mini- 
mum limits" sample. The corresponding estimates for the intercept 
of the metallicity-redshift relation are —0.65 ± 0.15 for the "max- 
imum limits" sample and —0.63 ± 0.21 for the "minimum lim- 
its" sample. Thus our data support the conclusions of Khare et al. 
(2004) and Kulkarni et al. (2005) that the global mean metallicity 
of DLAs shows at best a slow evolution with redshift. 



6 CONCLUSIONS AND FUTURE WORK 

The MMT observations presented here, together with our previous 
MMT and HST data have more than doubled the DLA Zn sample 
at z < 1.5 and more than tripled the sample at 2 < 1. Combining 
our data with data from the literature, we find that the systems with 
higher [Zn/H] also have stronger Cr depletion [Cr/Zn]. Analysis 
of the TVhi -weighted mean metallicity vs. redshift for our sample 
combined with previous data from the literature supports previous 
conclusions that the A^hi -weighted mean global DLA metallicity 
rises slowly at best and does not reach solar levels by z = 0. 

Questions still remain pertaining to dust depletion and selec- 
tion effects, and linked to this, the metallicity-redshift relationship 
of DLAs. Despite the large improvement from our surveys, DLA 
samples at z < 1.5 are still relatively small. Clearly, more observa- 
tions of DLAs at z < 1.5 are needed for improved statistics in this 
vast epoch (~ 9 Gyr). The large number of SDSS absorbers now 
available present a good opportunity for this purpose. 

Two of our systems have higher iVznii than the empirical up- 
per limit (log iVznii < 13.15) that has been attributed to dust ob- 
scuration in previous studies. The fairly high Zn abundances in 
these 2 DLAs suggests that metal-rich DLAs may be a rare class 
of objects, but can be found in a large enough sample. It would 
be interesting to confirm the metal-rich nature of these DLAs with 
higher resolution spectra in the future, and find similar other sys- 
tems by observing a larger sample. If a large number of such sys- 
tems are found, they may make a significant contribution to the 
cosmic budget of metals. 
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Figure 14. [Cr/Zn] vs. log ./Vjji for our data as well as measurements from 
the literature. Only detections are plotted. 



Figure 16. The metallicity [Zn/H] vs. z for the DLAs in this sample as well 
as from the literature. 
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Figure 15. [Fe/Zn] vs. log JVhi for our data as well as measurements from 
the literature. Only detections are plotted. 



DLAs may contribute as much as 17 % of the total metals. A(g - i) 
measurements are available from SDSS photometry for thousands 
of quasars with candidate DLAs and sub-DLAs. It would be in- 
teresting to obtain [Cr/Zn] and [Fe/Zn] for a large sample of such 
objects to study dust depletion patterns in quasar absorbers and to 
examine trends between quasar reddening and dust depletion. 
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It has been suggested that the low metallicity found in DLAs 
is caused by metallicity gradients. On the basis of a small sample 
of objects Chen et al. (2005) reported a difference between emis- 
sion line and absorption line metallicities. If metallicity gradients 
are present in DLAs, then metal poor DLAs may be more likely 
to be found because larger impact parameters are more probable. 
This may be because (a) they correspond to larger cross-sections 
and (b) they correspond to l ess dust. On the other ha nd, other stud- 
ies JSchulte-Ladbeck et all2005lJBowen et alj2005Tl have found no 
difference between emission and absorption line metallicities for 
other DLAs. Given that the metallicity gradients are not well estab- 
lished and fairly weak even in nearby spirals, it is not clear if metal- 
licity gradients fully explain the low metallicities seen in DLAs. 

So far, no DLA has been observed with a super-solar metal- 
licity. The key to finding such systems may lie in observing faint, 
reddened QSOs. Vladilo & Peroux (2005) have estimated that dusty 
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